Distribution of reflection eigenvalues in many-channel chaotic cavities with absorption 
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The reflection matrix R — S, with S being the scattering matrix, differs from the unit one, when absorption 
is finite. Using the random matrix approach, we calculate analytically the distribution function of its eigenvalues 
in the limit of a large number of propagating modes in the leads attached to a chaotic cavity. The obtained result 
is independent on the presence of time-reversal symmetry in the system, being valid at finite absorption and 
arbitrary openness of the system. The particular cases of perfectly and weakly open cavities are considered in 
detail. An application of our results to the problem of thermal emission from random media is briefly discussed. 
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Introduction.- When absorption in a cavity is finite, part of 
the incoming flux gets irreversibly lost in the walls, breaking 
thus the unitarity of the scattering matrix. The mismatch be- 
tween incoming and outgoing fluxes in the scattering process 
can be naturally described by means of the reflection matrix 
R = S}^Sj, where denotes the subunitary (at nonzero ab- 
sorption) scattering matrix to be precisely defined below. The 
matrix R has the positive eigenvalues Vc < 1, the so-called 
reflection eigenvalues. Their statistical properties in chaotic 
resonance scattering attract much attention presently, both ex- 
perimentally 1 1, 2] and theoretically A. ^ JL1..Sl.9\, using 
the random-matrix theory approach iTftlllll . 

In the limit of weak absorption, R was found to be re- 
lated to the Smith's time-delay matrix at zero absorption those 
statistical properties were extensively studied in recent time 
lll2lfl3lll4l[l5l fT^. Such an analysis was recently generalized 
by us to the case of arbitrary absorption |8], see also |9] for 
the related study. In the particular case of the single-channel 
cavity, obtained results explain partly the recent experiment 
101 on the reflection coefficient distribution. 

The opposite case of a large number of equivalent chan- 
nels plays an important role in chaotic scattering, correspond- 
ing to the semiclassical limit of matrix models ITtI \\m . In 
this case N strongly overlapping resonances (poles of the 5- 
matrix) are excited in the open cavity through the scattering 
channels, whose number M scales with N . A new energy 
scale appears then in addition to the mean level spacing A of 
the closed cavity 1 12, 18]: the empty gap between the real axis 
in the complex energy plane and the cloud of the resonances in 
its lower half plane. In the physically justified limit M /N^l 
llT2il (though both AI, N^oo) considered below this gap is 
given by the well-known Weisskopf width Tw = MT A/27r, 
with T < 1 being the transmission coefficient JItIi . In the 
presence of finite absorption the widths of the resonances ac- 
quire an additional contribution, the absorption width Fa 0]. 
The two competing processes, the escape into continuum and 
dissolution into the walls, are therefore characterized by the 
ratio Ta/Fw — l/T, where the convenient dimensionless 
parameter 7 = 27rro/AAf measures the absorption strength. 

In this paper we derive the probability distribution function 
'P{r) = A/^^ X^cli '^('' ^ "^c) of reflection eigenvalues, with 
the bar indicating the statistical average, in the limit of the 
large number of channels. This function has, in particular, 
an important application for the statistics of thermal emission 



from random media and is known only at perfect coupling 
(r=l) |5]. Here, we calculate it at arbitrary T. 

Our starting point is the general relation established in 
Ref. 1 8] between R and the effective non-Hermitian Hamil- 
tonian H = H — {i/2)VV'' of the open system, with the 
Hermitian part H standing for the closed counterpart and the 
amplitudes describing the coupling between N intrinsic 
and M channel states. It is achieved by making use of the 
equivalence 0,0,01 of uniform absorption to the pure imag- 
inary shift = E + ^Ta of the scattering energy E, thus 
giving = SiE-^) = 1 - iV^E^ -Uy^V. This leads to 
the following connection between the reflection matrix R and 
the time-delay matrix with absorption Q-^ ifi — 1) 01: 



i? = 1 - TaQ. 



(1) 



where = Q{E^) = _7^)t]-i(£;^ The 

distribution function 'P(r) is therefore directly related as 



p(r)=7-ip, (7-1(1 -r)) 



(2) 



to the distribution ^^(r) = M'^ YlcLi ^('^ ~ ^^cltu) of 
the proper delay times (eigenvalues qc of Q^) scaled with M 
and measured in units of the Heisenberg time tn = 27r/A. 

The function (|2} is normalized to unity. The average reflec- 
tion coefficient (r) = M^^tr i? is given by 



= / dr rV{r) = 
Jo 



r + 7(i-r) 

T + 7 



(3) 



This result follows from the connection (Sl between (r) and 
the "norm-leakage" decay function P{t) introduced in |T^: 
(r) = 1 - 7[1 - Ta dte-^-*P{t)]. In the limit considered 
P{t) reduces to the simple exponent e^^^" *, since the widths 
ceased to fluctuate iflMllQI] . Below we will also derive Eq. Q 
using a different method. 

The saddle-point equation.- The scaUng limit studied is far 
from being trivial and can not be obtained from the known ex- 
pressions for at finite M |81 simply by letting Af— >oo 
there. One needs to start from the original definition that de- 
termines the seeking distribution through the jump of the re- 
solvent G{t) = M-hr {t-TaQj)-^ = l/t+ {-f/t^)K{t) 
along the discontinuity line t > as follows: 



V(r) 
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7r(l 



;lmK{t-iO)\^^^_ 



(4) 
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Following Refs. I81 II6II . one can obtain the exact (at iV— >oo) 
representation for K, which is convenient for our purposes, 



K{t) = {^StT 



A - ii^=i^Ai ) ka 



(5) 



with the shorthand ((. . .))c / d^i{a){. . .)e(*^/2)^ and 

C = {j/t) str[(l - t/2)Aa ~ Wl - tAia] 

-strln[l + (2g)"i(A(T + (TA)] , (6) 

as the integral over the noncompact saddle-point manifold of 
the 8x8 supermatrix a subject to the constraint &^ = 1. A, Ai 
appearing above are the supermatrix analog of the Pauli ma- 
trices CT3,f7i and k = +1(— 1) in the space of commuting 
(anticommuting) variables. Definitions of the superalgebra as 
well as explicit parametrization, which depends on whether 
TRS is preserved or broken, can be found in flolEoll . At last, 
the constant g — 2/T— l>lis related to the transmission 
coefficient T. 

In the limit M^oo the integration over a can also be done 
in the saddle-point approximation. Function ^ is given by a 
saddle-point value of a, which is found, as usual, by equating 
the variance 6£ to zero. The structure of the latter equation is 
determined by that of C. A careful analysis shows that, inde- 
pendently on TRS, the supermatrix a reduces (and so does the 
corresponding algebra) in the saddle-point to the 2x2 usual 
matrix cr = aa^+bai with the imposed constraint +6^ = 1. 
We find that K — a — — </2)/ ^/l — t, whereas the equa- 
tion SC = gives {-//t)VT^K + ib/{g + a) = 0. EHmi- 
nating in this equation a and 6, one arrives finally after some 
algebra at the following saddle-point equation 



1 



1 1 



iK + g)/2 



t 2 jK ±^[K + 5)2 + 1 + 4/72 - g2 _ 2/7 



(7) 



where the sign "+" must be taken. Independently on the taken 
choice of the sign, equation Q can be further equivalently 
represented as the following forth order equation in K: 



K-^ + 2gK^ 



(1 



+ -r 5 + - ' 

4 7 
K , 2, 



Sg 



(1 



1 

- = 0. 

T 



(8) 



In the zero absorption limit, 7^0 with fixed t/^ = t, this 
equation simplifies to the cubic one found earlier 1 16]. 

The choice of the sign "+" is justified by two reasons. First, 
only then the complex solution of Eq. Q yields the distribu- 
tion 7'(r) nonzero at positive r. Second, the function K{t), 
as follows directly from its definition, must have at t— >— cx) 
the positive limit Koo ~ i//^tr < 1 |8], which is the mean 
scaled proper delay time (r). The equation for the latter fol- 
lows then readily from (jSjl as [{g + 2 /j)'^ - 1] (r)2 + 4(t) - 
4/72 = that gives (t) = T/{j + T), reproducing thus 
Eq. (|3 by virtue of (r) — 1 — 7(r). 

The higher moments of the distribution can be easily calcu- 
lated in the same way, exploiting the 1 /i-expansion further 
For example, substituting K — (t) + ^{t^) /t + 0{t^^) in 



(|8} and equating there the 1/^-term to zero, we find (t2) = 
(r)2[72 + 2(7 + r)]/(7 + T)2. That yields the variance of 
reflection eigenvalues as follows: 



72r2(27 + 2r(i^7)-r2) 
(7 + r)4 



(9) 



Remarkably, the variance and, therefore, fluctuations of r are 
suppressed in the both limits of weak and strong absorption. 

To understand qualitatively the structure of the solutions 
of the saddle-point equation, it is instructive to consider the 
behavior of 1/t as a function of K shown schematically on 
Fig. 1. We readily see that the distribution function 7^(r) is 
nonzero only in the finite domain rmin < r < Tmax- The 
value of the upper border Tmax < 1 at any g and 7. In a vicin- 
ity of Tinax the distribution behaves as P(r) cx V^max — r. 
For the lower border we find rmin — and Vir) cx l/v^, 
when (5 — 2/7)2 < lor absorption being from the interval 
T < 7 < T/Yl - T), and r^in > with P(r) cx ^ - ^nin 
otherwise 12111 . Further analytical study is possible in the fol- 
lowing particular cases, which we consider in detail now. 

Perfect coupling, g = T = 1. In this case, the region be- 
tween the two vertical dash-dotted lines on Fig. 1 shrinks to 
the line K ^ ~1. Eliminating the resulting common factor 
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FIG. 1: The schematic plot of 1/f as a function of K. The solid 
and dashed lines correspond to the choices of and "— " sign in 
Eq. 0, respectively. Altogether they form the complete set of so- 
lutions of the forth-order equation js). In the shadowed region the 
latter equation has two complex-conjugated roots (and two real roots 
irrelevant for our purposes), with t^^ > 1 (r > 0). Thus, only in 
this region the corresponding density is nonzero. Another solution 
with < < 1 is unphysical. 
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{K + 1), one gets from Eq. (|8} the following cubic equation 



1 - r 



-1-7 



= 0, (10) 



written already in the variable r = 1 — t. We find that the 
discriminant D = [(l+7)/(37'^r)]'^(r+ — r)(r — r_)(l — r)"' 
of this equation, with 



r± = 



8 + 2O72 -74±7(8 + 72)3/2 



(1 + 7)3 



(11) 



is positive at ?■ > (and thereby ImK is nonzero) only in the 
domain max[0, r_] < r < r+. This sets explicitly the exact 
borders for the searched distribution. 'P{r) can be found from 
ilOi by applying Cardan's formulae, reproducing exactly the 
result of Ref. |5] obtained by a different method. When ab- 
sorption is weak, 7 < 1, the behavior of 7'(r) can be well ap- 
proximated by the following simple interpolation expression 



27r (l-r)2^ 



(12) 



with = 1 



+ ^7 + 0(7 ) being the normalization constant. 
In the opposite case of strong absorption, the distribution is 
found to be close to 



7>1 



(r) 



(1 



(13) 



The limiting distribution J12> or (I13t becomes asymptotically 
exact as 7 diminishes or grows, respectively. All these fea- 
tures are illustrated on Fig. 2a. 

Special coupling, — 1 + 4/7^. Under this condition 
Eq. reduces to a quadratic equation in K, giving readily 



2v/r^v/rg/r-l 
TTrg (1 — r)2 



(14) 



with r, = 2/(1 + V1 + 7V4) = 2/(1 + g/vV^)- 

Weak coupling, 5 « 2/T ^ 1. One can find approximate 
expressions for Vir), making use of the scaling K^O{\/ g) 




FIG. 2: The exact reflection distribution (solid lines) compared to 
the corresponding approximate expressions <12t . <13> . and <15> at 
perfect (T=l, left) and non-perfect (r=0.5, right) coupling. 



in the exact equations. When absorption is weak we may ex- 
pand the square root in and keep only the main contribu- 
tion (2g/^){l + 2gK)^^ instead of the last term there. That 
leads again to a quadratic equation in K, giving readily the 
behavior 7'(r) w /^gn^y/^ ^r^^^^ — rj (1 — r)^ near the 
upper border Tmax ~ 1 — 7/83. For 7 from the small inter- 
val [T, T/(l - T)], when 7min==0, V(r) is given in the lead- 
ing approximation by Eq. WW . with replaced by rmax, be- 
cause the applicability condition for WW is now satisfied up 
to 0{g^^). At other values of small 7 there is no reliable ap- 
proximation for the lower border available and rmin > is 
to be found from the general equation 1 2 1 ] . At last, for very 
strong absorption and arbitrary non-perfect coupling, T ^ 1, 
we can use scaling K^0{1/^) in 0. That allows us to 
replace there the last term with ^{K + g)/{l + gK — 2/7), 
yielding finally 



7»1 



(r) = 



C v/(nrrAf^^7)(7^^7^TA^ 



27r (l-r)2(l-T + r) 
with the normalization constant C«7 — 2(1 — T) and 



(15) 



1 -T 



^lii^ and A..2V2rji^ 



7 



7 



that is valid when ro ^ Ar or 7 KT"^ / (1 — T). Expression 
(I15> becomes asymptotically exact as 7 grows, approaching 
the "semi-circle" distribution with the center at and the ra- 
dius Ar, see Fig. 2b. 

Thermal emission.- As an application of our results we 
consider thermal emission from random media. In his sem- 
inal paper Beenakker |5] has shown that the quantum opti- 
cal problem of the photon statistics can be reduced to a com- 
putation of the 5-matrix of the classical wave equation. In 
particular, chaotic radiation may be characterized by the ef- 
fective number v^ff degrees of freedom as follows: v^s/v = 
(1 - (r))V((l - r)2) < 1 15], with v^n = v for blackbody 
radiation ll23ll . We find using (|9} that at any 7 and T 



V 



il + Tf 



(t2) -^2 + 2(^ + T) 



(16) 



and a mean photocount n = vf^/T/ (7 + T), with / being a 
Bose-Einstein function. The earlier result |5] is reproduced 
at T—1. Upon the substitution of 7 with —7 one can use 
[5] n and v^f^/v ( I16> even for amplified spontaneous emis- 
sion below the laser threshold, 7 < T. In this case our mean 
photocount agrees with the findings of Ref. 0|, where the 
general theory of photocount statistics in random amplifying 
media was developed. In the limit of vanishing absorption or 
amplification, the ratio v^a/v is T/2. The large 7 expansion 
1 - 2(1 - r)/7 + 0(7"^) shows that the saturation 
to the blackbody limit gets slower when transmission T < 1. 

In conclusion, for many-channel chaotic systems we have 
derived the general distribution of reflection eigenvalues at ar- 
bitrary values of absorption and transmission. We note that 
due to a duality relation |5, 25], an amplifying system in the 
linear regime (Fa < ^w) is directly linked to the dual ab- 
sorbing one through the change of the sign of Fa in Q and 
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correspondingly thereafter As a result, the reflection matrices 
(and their eigenvalues) of dual systems are each other's recip- 
rocal. Therefore, the analysis presented can straightforwardly 
be extended to the case of linear amplification that might be 
also relevant for the rapidly developing field of random lasers 
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